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Abstract 

We compute the two-loop virtual corrections to the flavour changing neutral cur- 
rent process b — > dt + 1~ . As calculation techniques integration by parts identities 
and the method of differential equations are used. The result is presented in closed 
form as a function of q 2 /m 2 , where q 2 is the invariant mass of the lepton pair and 
nib is the 6-quark mass. 



1 Introduction 



Flavour changing radiative B-decays play a crucial role in testing the standard model, 
because they proceed entirely through loops and are therefore sensitive to the propaga- 
tion of virtual new particles. At present, one of the most stringent constraints on new 
physics from rare S-decays arise from the inclusive decays B — > X s 7 and B — > X s £ + £~ 
which have been studied in great detail in the literature [U 13 El • Within the present 
uncertainties, the experimental results jUEHO] f° r both decays are fully compatible with 
the theoretical predictions IHj arising in the SM. 

Concerning the two-body decay B — > X,7 calculations are complete at NLO in 
QCD 7j. For the process B — > X s £ + £~ semi-numerical results are available for the low 
and high range of q 2 (invariant mass of the lepton pair), i.e. after cutting out the nonper- 
turbative uu and cc resonances [Hj. There only the numerically small Wilson coefficients 
arising from penguin contributions were partly neglected. In these two decay modes, 
the branching ratio and the invariant mass distribution/forward-backward-asymmetry 
for the latter one are important observables that can be used to put constraints on 
parameters of models beyond the SM. 

The transitions b — > 57 and b — > s£ + £~ have, to a good approximation, vanishing 
CP-asymmetry, because the CKM-element V* s V u b can safely be neglected and therefore 
the CP-violating phase drops out of all observables. This changes when one deals with 
the b — > d transition, because V* d V u b is of the same order than the other two matrix 
elements contributing to the process. Therefore we have, in addition to the before 
mentioned observables another one, the CP-asymmetry which can be used to test the 
SM or to constrain the angles of the CKM-triangle. 

As a consequence of the new CKM-structure, at NLO one has to compute in addition 
to the two-loop diagrams with an internal c-quark propergating j^J, the ones where the 
c-quark is replaced by a w-quark. We present the non-factorizable two-loop virtual 
corrections in an analytic form for arbitrary q 2 , where q 2 is the invariant mass of the 
lepton pair. 

For the decay mode with a real photon in the final state, the results can completely 
be extracted from [7j. In that paper the diagrams were calculated for a whole range of 
operators in and beyond the standard model. Here we restrict ourselves to SM operators 
and furthermore we neglect the penguin Wilson coefficients because of their smallness. 
Recently Asatrian et.al. ^H] presented the results for B — > X < ji + t~ where the same non- 
factorizable two-loop corrections were calculated. They have presented it as a power 
series in q 2 /m 2 , as we give analytic formulae for the corresponding diagrams. After 
expanding the results given in this paper, we fully agree with their findings. 

The phenomenological analysis of the decay will be done in jTHi. so we restrict our- 
selves here to the pure computation of the two-loop diagrams. 

The paper is organized as follows. In the next section the effective Hamiltonian for 
the considered process is given. Section |H] contains the description of the method used 
to reduce the two-loop integrals to a few scalar Master Integrals (Mis). The following 
section deals with the calculation of these integrals. Then we give the results for the 
virtual two-loop corrections and finally some concluding remarks are given. In the 
appendices the results for the Mis, the unrenormalized virtual corrections and the one- 
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loop integrals needed for renormalization are given. 



2 The Effective Hamiltonian 

In the Standard Model the weak effective Hamiltonian forb—>-d transitions is given by 
AG F 



V2 
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,q=u,c i=3 



+ h.c, (1) 



where X g d) = V* d V qb . 

In order to match our computation with the NNLL results for the Cj, we use the 
operator basis 



0\ = (d L7fl T a qL )(q L YT a b L ) } 0\ = (d Llfl q L )(q L Yb L ), 

3 = (d L 7m b L ) (Q^Q), 4 = (d L 7 M T a b L ) Y,, (QY T a q), 

O, = (d Ll ^ p b L )J2 q (qrYYq), O e = {d L ^ lvlp T a b L )J2 q (Ql*ri p T a q), 

7 = (d L a^b R ) F, u , 8 = ^ (d L a^T% R ) G a , 

9 S 9s 

9 = ^f (d L 7 „ b L ) Y,e (^ t), O w = (dL 7m b L ) (^75 i), (2) 

9s 9 s 

where the sum over q and i extends over all light quark and lepton fields, respectively. 
gem(gs) is the electromagnetic (strong) coupling constant, qi and qn are the left and 
right chiral quark fields, F^iG^) is the electromagnetic (gluonic) field strength tensor 
and T a are the color matrices. 

The one- and two-loop matching conditions for the Wilson coefficients at the scale 

CM = CfW) + ^ cfW) + (^) 2 Cf W) + 0(a* s ) (3) 
can be found in whereas the anomalous dimension matrix 

^^M£)>+(£)> ,+ °w> < 4 > 

which is needed to evaluate the coefficients down to the scale /i ~ nib has been calculated 
to NNLL order in [IE] . 

It has been shown that decay amplitudes for the decay of a 5-meson into light 
mesons in the kinematic region of large recoil of the light meson can be systematically 
computed in terms of form factors, light-cone distribution amplitudes and perturbative 
hard scattering kernels jTTl HH1 E2|- Schematically, the amplitude can be represented 

as mi 

(£ + £-p a \H eS \B) =C a £ a + $ B ®T a ® $ p , (5) 
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Figure 1: Two-loop contributions to {y*d\H e s\b) that have been calculated in this paper. 
The circled crosses mark the possible insertations of the virtual photon line. 

where a =_L, || refers to a transversely and longitudinally polarized p-meson, respectively. 
In this equation £ a represent universal heavy-to-light form factors ITU] and $ light- 
cone-distribution amplitudes. The factors C a and T a are calculable in renormalization- 
group improved perturbation theory. In this paper we compute the NLO non-factorizable 
two-loop contributions to the C a . The corresponding diagrams are shown in fig. ^ In 
addition, there are factorizable two-loop vertex corrections. As the difficult part of 
the calculation lies in the non-factorizable diagrams , we refrain from presenting the 
factorizable ones. For the exclusive decay these have been calculated in [TJ] whereas we 
refer to JU] concerning the inclusive mode. 

3 Reduction to Master Integrals 

In this section we present the calculation of the two-loop matrix elements {d^*\0\ 2 \b), 
were 7* denotes a virtual photon. From these we easily get the quantities (d£ + £~\Oi 2 \b) 
we finally wish to compute. The corresponding diagrams are shown in figffl 

It is always possible to decompose the matrix element into the following form: 

(j*(q,n)d(p')\H eS \b(p)) = (d(p')\dX»b\b(p)} = F^(q 2 )(d(p')\d(l +1 5 )q»b\b(p)} 

+F^(q 2 )(d(p')\d(l + 1 5 )a^q u b\b(p)) 

+F^(q 2 )(d(p')\d(l + 7 5 )(gV " ff{)bW)) (6) 
The scalar form factors F^ q \ F^ 7 \ F^ are obtained by taking the trace 

F«(g 2 ) = Tr(iTX M ) (7) 

with the projectors 

Pt = {j + m b ) (C tl <f + C l2 p" + C l3 7") f (8) 
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and some scalars Cy. The factors (tf + mj) and f account for the onshell condition 
of the b— and s— quark lines, respectively. can give a nonzero contribution to 

individual diagrams but vanishes in the sum of the diagrams because of electromagnetic 
gauge invariance. This gives a good check of the computation which is described in the 
following. After projecting out the form factors, these are of the form 

where the following scalar products and propagators containing loop momenta appear: 

Si = k 2 , iS>2 = I 2 , S3 = k ■ I, Si = k ■ p, 

S 5 = k-q, S 6 = l-p, S 7 = l-q, (10) 



T>, = — Vo = V-i = — 

k 2 2 {k + p-q) 2 ' 3 {l + qf 

Va = 4, Vr = ~r, r- Va 



I 2 ' {l + k + q) 2 ' (k - p + q) 2 - m 2 ' 

Vl = (l + k) 2 ' Vs = {k + p) 2 ' Vd = {k-p) 2 -m 2 ' ^ 

For the diagrams in fig. ^ we get at most five powers of the propagators V% and three 
powers of the scalar products S^ 

The next step of the calculation is to reduce the above integrals (in our case there are 
a few hundred) to a few Mis by means of the so-called Laporta algorithm |23] 124"] using 
Integration by Part (IBP) [213 E] and Lorentz invariance [22] identities. IBP identities 
use the fact that in dimensional regularization we have the following equalities: 



ddhd d l— V" 12 ' 7 \ ~ 



.:•.%■• ) =° (12) 

where t> M = g M . For the kinematics in the present process the Lorentz invari- 

ance identity for an arbitrary scalar integral takes the form 



d d kd d l 



u 9 .. d 2 ,. d 2 .. d 

P'lP TTH-P-II TTH + I v-z—-p q^jr- b 
ap^ aq^ aq^ ap^ 



QU\ QT12 QUJ 

Q 1 Q 2 . . . o 7 

1?mi<r>ni2 Ti m 9 
' 1 '2 • • • ' 9 



0.(13) 



This way we get nine identities for each scalar integral. These identities will contain 
other scalar integrals with different numbers of scalar products and propagators. After 
making up the identities for all the appearing integrals, the aim is to solve these linear 
equations for integrals of a simpler type. In principle this is a simple task but with the 
diagrams considered here we get linear equation systems up to about 3200 equations and 
therefore some organization is required to solve them. 
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a) b) c) d) 




Figure 2: Scalar Master Integrals appearing in our calculations. The dashed lines denote 
massless propagators whereas the solid lines represent the ones with mass m&. Dotted 
propagators have to be squared. Wavy/solid/dashed external lines are onshell by the 
amount g 2 /m^/0, respectively. 



As stated above, we start with integrals with at most five powers of different propa- 
gators. It is possible to express all of these in terms of integrals containing four different 
propagators. Concerning these simpler topologies, there exist some integrals where no 
equation relating these to simpler topologies can be found by that method. The same 
happens for integrals with three different propagators. The remaining integrals are the 
Mis which we solve by the method of differential equations. Fig. El lists all the Mis 
appearing. 

4 Calculation of the Master Integrals 

The last step in the calculation of the form factors is the evaluation of the Mis. We find 
five Mis with three propagators and seven with four different propagators, where in the 
latter case, two of the propagators show up squared. The 3-denominator topologies in 
fig. Eh-1 and the MI in fig. |2k which factorizes into two one-loop integrals can be calcu- 
lated trivially. Nevertheless we have given the results for completeness in appendix [XJ 
like the results for the other Mis. For the remaining Mis we use the method of differential 
equations (251 ESI I2H| ■ In our case, the equations take the form 

^me,mi) = ^*(£-£) Mfe 2 x> a*) 

where the derivative is understood to act on the integrand of the corresponding MI. After 
evaluating the right hand side of (jl4)L there will be terms containing the MI itself but 
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also terms with other scalar integrals. We now use the same method as described in the 
last section to express these in terms of Mis. We end up with two different possibilities. 

In the first case the derivative with respect to q 2 of the MI gives terms proportional 
to the MI itself and some Mis of a simpler topology, where one propagator is completely 
canceled, not just reduced by one power (figs. EJ>c). In this case we are left with an 
inhomogeneous linear differential equation with 3-denominator topologies as inhomo- 
geneities that we have already solved. Solving this equation is even easier by expanding 
the Mis in e: 

CO 

M\{q\ml)= ^-MI% 2 ,m 2 ) (15) 

t=-2 

After expanding the differential equation (the coefficients are e-dependant) in e, we 
can solve it order by order in e, the inhomogeneous parts either coming from the 3- 
denominator topologies or from the part of the 4-denominator MI in the e-expansion 
that has already been solved. 

The second possibility arises when solving the Mis that are shown in figs. Oi-e and 
figs. |2F-g, respectively. After reducing ()14j) to contain only Mis, we are left with a coupled 
inhomogeneous differential equation for the two 4-denominator topologies. When we 
expand this system in s, the equations decouple order by order. In both cases we used 
Eulers variation of the constant to solve the different parts. 

This method can be used in principle to get the Mis to any desired order in e. As 
described in the previous section, the form factors are expressed in terms of Mis with 
certain coefficients, which can depend on e. Since we want to know them up to the finite 
parts in the limit e — > 0, it turns out that we do not have to exceed the second order in 
the expansion for the Mis. For one of the integrals (fig|2t) we actually only need to go 
to the 1/e-term. 



5 Renormalized matrix elements 

In this section the main result of this paper, the MS-renormalized NLO non-factorizable 
two-loop virtual correction relevant for the decay b — ► d£ + £~, is given. We present the 
matrix-elements as in [TOj, that is after the decomposition 

(den- \o?\b) non _ {act . = (F^(d 7 ) tiee + F^(d 9 ) tree ) (i6) 

two-loop Q-T^ ^ ' 

where (0 7 ) tree and (Og) tTee denote the tree level matrix elements of 7 and Og, respec- 
tively. Furthermore we have defined the rescaled operators 

^7/9 = ^07/9. (17) 

Note that our sign in the above definition differs from that in |lUj . 

As we use the MS-scheme renormalization constants, all unrenormalized 1-loop form 
factors are assumed to be multiplied by {An)~ le e lei . They are given in appendix IHl and IO 
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Since the Wilson coefficients C^-q for the penguin contributions were neglected, the part 
of the renormalized effective Hamiltonian which contains the counterterms needed for 
renormalization reads 



i=l 



2 10 12 

E + E 5Z v°i + E <^A W 

J=l 3=3 j=ll 



+ .... (18) 



The operators Oi-io are given in (j2j), whereas On/12 are evanescent operators vanishing 
in d — 4. In order to match our computation with the Wilson coefficients of [16J, we 
choose them as 

On = (d Ll , lulp T a u L )(u Ll ^YT a b L ) - 160?, 
^2 = K7m7,7p^)(^7 m 7V&l) - IQOl (19) 
The operator renormalization matrix Z — 1 + <5Z is given by ^5] 



47T 



a. 



(X q I e 3s 9e 27s 12e 9e 



OOfOOOO-^Oi 

3e 9e e 







58 n 64 , 1168 



4rr 



243e 729e 1 243e 2 



+ S +0(af) (20) 



000000 ^ o + f + 148 



81e 243e 1 81e 




As a lot of entries in this matrix are zero, we only have to compute the one-loop matrix 
elements of 0\ 2 4 n 12 an d furthermore the tree-level matrix elements of are needed. 
Mass and wave function renormalization is a higher order effect and we do not have to 
include it here. The strong coupling constant in the definition of the operator Og has to 
be renormalized where Z 9s is given by 

z _ Osfol _ 2 

9s 4tt 2 e' 3 ; V ; 

Note that there is no renormalization needed for a s in O-j because mixing of the four- 
quark-operators into Oj vanishes at one-loop and therefore the corresponding compo- 
nent in 5Z only starts at 0(a%) and the coupling constant renormalization for O7 is a 
higher order effect. 

As only two different color structures arise in the different diagrams, we can split the 
form factors up into 

*fi = M*), (22) 
ifj = -6A(s), (23) 

F^ = B(s)+4C(s) } (24) 

F^ = -6B(s) + 3C(s), (25) 
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where the functions A(s), B(s) and C(s) are given below. 
The following definitions are used in the formulae: 



g 2 , 



4m\ 



m 



2 ' 



(26) 



1 * / T 



X 3 



+ 



1 1 I T 

x 2 = ---v / ^l, 



X4 



2 2 V / ^T' 2 2 V / ^T' 

/i ~ m b denotes the renormalization scale, ( the Riemannian Zeta function and 

ln(l -t) 



(27) 



U 2 (x) 



dt- 



ps) 



is the Dilogarithm. 

The functions A(s), B(s),C(s) are as follows: 



, . , 104 .ml. 



4s 



243 > 2 27(1 -s) 
1 



729(1 - s) 2 
2 

243(1 - s) 3 

+18iris(l - 2s) ln(s) 



Li 2 (s) + ln(s) ln(l-s) 
6s (^29 - 47s) ln(s) + 785 - 1600s + 833s 2 + 6ni (20 - 49s + 47s 2 ) 



2V-2 - l( - 4 + 9s - 15s 2 + 4s 3 ) arccot(v / ^ - 1) + 9s 3 ln 2 (s) 



2s 



243(1 - sY 



36 arccot 2 (v / ^- 1) + tt 2 - 4 + 9s - 9s 2 + 3s 



(29) 



B(s) 



22 2 
(4 - 34s - 17vns) ln(^) + 8sln 2 (^) + 17sln(s) ln(^ " 



243s 
(2 + s)v / ^I 



729s 



2 

- 48 ln(^) arccot(v / ^ - 1) - 18vr ln(z - 1) + 3i ln 2 (z - 1) 



-24iU 2 (-x 2 /x l ) -57r 2 i + 6^-91n 2 (a;i) + ln 2 (x2) -21n 2 (x 4 ) 

+6 ln(xi) ln(x 2 ) — 4 ln(xi) ln(x 3 ) + 8 ln(xi) ln^)) 
-12?r (2 ln(xi) + ln(x 3 ) + ln(x 4 )' 



243s(l - s) 



4s ( - 8 + 17s) (Li 2 (s) + ln(s) ln(l - s)) 



+3(2 + sjfo - s) ln 2 {x 2 /x 1 ) + 12vr( - 6 - 5 + s 2 ^arccot(V^ = 1) 
- 18s (l20 - 2115 + 73s 2 ) ln(5) 



21875(1 - s) 2 

-288 -81 + 934s 2 - 092^ ! ISnis { 82 - I 73* 73s 2 
4 



2V 7 ! 3 ! (4 - 35 - 18s 2 + 16s 3 - 5s 4 )arccot(v / ? 3 T) 



243s(l - s) 3 

-9s 3 ln 2 (s) + 27ris ^8 - 33s + 51s 2 - 17s 3 ) ln(s) 

72^3 - 8s + 2s 2 )arccot 2 (v / ^ - 1) 



7295(1 - s) 4 



vr 2 54 - 53s - 286s 2 + 612s 3 - 446s 4 + 113s' 



m s 16 , / s s 428 64 16 . 

C(s) = + CCS) + — 1T1 

v ' 81 V 243 27 SV ' 81 



(30) 



(31) 



In fig. El we give a comparison between the numerical values of the NLO corrections 
to the form factors , F± , F2 as computed in this paper and the result of the ex- 
pansion to third order in q 2 /ml which agrees with the result presented in ^U]. The 
physical "threshold" of the charmonium resonances begins at M 2 JfiJ w 9.6 GeV 2 , but the 
perturbative approximation is expected to fail earlier. Model-dependent studies of the 
charmonium contributions suggest that the perturbative approximation should be valid 
up to q 2 ~ (6 — 7) GeV 2 |2Z1 EHj- As can be seen in fig. 01 the approximation up to 
that value of q 2 is fairly good but fails as expected in the high q 2 region, beyond the 
resonance. The largest deviation (at q 2 = 7 GeV 2 ) comes from Im(F^) and amounts to 
1.3%. 



6 Concluding discussion 

In this paper we have presented the non-factorizable two-loop virtual corrections con- 
tributing to the process b —> d£ + £~. We used the technique of integration by parts for 
the reduction of the few hundred integrals to a few Master Integrals. To solve these Mis 
we apply the method of differential equations which allowed us to calculate the two-loop 
diagrams analytically. Asatrian et. al. ^U] presented these calculations as an expansion 
in q 2 /m 2 , whereas our result is valid for the whole range of q 2 . When expanded in q 2 /ml, 
our result agrees with the calculation made by Asatrian et. al. HT 
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Figure 3: Real and imaginary parts of the form factors F± , , F 2 . The solid line 
represents the exact result to 0(a s ), whereas the dashed line shows the result expanded 
in q 2 jm\ up to the third order. 



A Solutions for the Master Integrals 

In this appendix we give the solutions for the three- and four-point Mis needed for the 
calculation of the form factors for the two-loop integrals. Especially the three-point 
integrals and the four-point integral shown in fig. |2k are easy to calculate but are listed 
here for completeness. 
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In the following the integrals are normalized according to |29j : 

^ - wt-m (32) 

We give the results up to the required order needed for our evaluation of the form factors. 
A.l 3— denominator topologies 



+ 7 - |vr 2 - 3 ln(s) + - ln 2 (s) + 3vri - iri ln(s) + 0(e) \ (33) 



'WWV -T\/\/\/\/ 



fi 2 ^ I {d d k}{d d l} ' 



V{P & V 



7 



+ |g ^2595 - 364vr 2 + 1380ttz - 48vr 3 i - 240C(3)) + 0(e 3 ) J 



/i 2 ^ / {d d A;}{^/} ' 



V 2 V 6 V 9 

2\ -2e 



P 2 P 5 P 6 



2 /m 2 \ 2E ( 1 3 7 1 2 ■: 
2^ + 2i + 2 + r 



(34) 



+ |(l5 + 7r 2 -2C(3))+0(, 2 )} (36) 



V 2 V 6 V 8 

vwv\| rww (s — ► (2^) 
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A. 2 4— denominator topologies 



^ J {d d k}{dH}^ 



2 

/i 2 ^ \e 2 e v ~ ■••-/■ g 



s \ 2e f 1 2 , 7 



+ "(2 + 7TZ) + 12 - -7T 2 + 8ttz 



2 ( 16 - y 7T 2 - 2C(3) + 12m - 7i 3 i ) + 0(e 2 ' 



(38) 



+ % ( 195 + 57r" N 



^V 2£ I ^ + 1 + - + - 2 + , 3 ,fln 2 (g) - 2mHs 
V/i 2 y \2e 2 2e 6 6 2(1 - s) \ KJ V ' 

12sLi 3 (s) - 12sYln(s) - 7r^Li 2 (s) 



6 V J 6(1 -s) 

- s ln(s) (W 2 + 15 ln(s) - 5 ln 2 (s) + 6 ln(s) ln(l - s)) 

+ 3ms ln(i) (4 ln(l - s) + 10 - 5 ln(i)) - 2svr 3 i + 6 (l - 3s) C(3) + 0(e 2 ) J 



(39) 



where 



a; 



Li,(.r): / (40) 



o 



is the Trilogarithm. 



Ai 2 ^ / {d d k}{d d l} 1 

-(?r{ai^(-^ + «) 



1 19 4 

+ Li 2 (s) - 5 ln(s) + - ln 2 (s) + ln(s) ln(l - s) + tt 2 - 2m ln(s) + 5vu + 0(e) 

2 2 3 



(41) 
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/r- i {d d k}{d d l} 



L 2(A-d) 
,2\ - 2e 



1 



v 2 v 6 v 7 v 9 



+ 



6(1 - s) 



/ m 2\ - ze r i i / \ 19 

= - (-*) 1 2? + 27 ( 5 " 2 ln(5) + 2?r V " 5 ln(5) + y - ^ 2 + 5m 

6Li 3 (s) -6^1n(i) - vr^Li 2 (s) - 3 ln 2 (s) ( ln(l - s) + 21 - 1 
+ ln(s) f 7T 2 + 6m ln(l - s) + 12ms - 6m) - n 3 i - 6C(3)1 + 0(e) J> (42) 



V 2 (P 6 ) 2 V 7 V 9 



1 /m 



2\ -2e 



ln(s)( ln(s) -2vri 



2(l-s)e~* v ' / V~ v '' ~"V 6(1 -s) 
18(ln(s) - vr^Li 2 (s) + ln(s) (5 ln 2 (s) - 91n(s)ln(l - s) - 7tt 2 

- 15tt2 ln(s) + 18vri ln(l - s)) - 3vr 3 i - 18C(3)1 + 0(e 



18 Li, (si 




2(4-d) 



v 2 v 5 v 6 v 9 



= ~{f) 2£ {^-^(4v^larccot(v^T)-5 
;f24Li 2 (x 2 /x!) - 12Ii 2 (-x 2 /x 1 ) - 15 ln 2 ^/^) + 121n 2 (x 2 ) 



12 



12 ln(xi) ln(x 2 ) + 24 ln(xi) ln(x 4 ) - 24 \n(x 2 ) ln(x 4 ) - 5tT 



114 

+ TTT + 



+ 6vr ( ln(xi) + ln(x 2 ) - 2 ln(z - 1) ) + 120 arccot(v / 2 - 1) 
1 



12 12(1-5) 



3-s) (W^M) 
- 36tt arccot( Vz - 1)) + tt 2 ( 33 - 13s ) + 0{e) 



(43) 



(44) 
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— — - ( 36 arccot^v 7 ^!) - vr 2 ) + 0(e Q 
1 — s)s V / 



(45) 



6a(l - s) 

B Unrenormalized matrix elements 2 |^) 

We present the results of the matrix-elements as in section 03 After splitting up into the 
different diagrams and the different color structures that can arise, one gets 



*l,u 



7 %-(7) 

2,11 



2\ -2e 



*5 



%-(9) 



2\ -2e 



7 %-(9) 

2,M 



^ 2£ (i? fis -( S )+4^( S 

■65^(5) + 3C fig -(s)), 



2\ -2e 



(46) 

(47) 
(48) 
(49) 



where the functions A g '(s), S g '(s) and C g '(s) for the individual group of diagrams are 
given in the subsequent paragraphs up to the zeroth order in e. 



Diagrams shown in fig. [Ql 



27e 9 27 
2s 



TYl 



6 Li 2 (s) - 3 ln 2 (s) + 6 ln(s) 1 + ln(l - s) - tt 



81(1 - S 

4 2 / \ 1 20 . 16 ., _ 

' = 27^ " We l 4 ln(5) " 5 " Am ) + 9 + 27 m ~ 27 ™ Ms) 



(50) 



81 



^— ^ [l2(l - 2s) (hi 2 (s) + ln(l - s) ln(s 

+6 ln(s) f hi(s) - 9 + 5s) - tt 2 f 15 - 171 



(51) 
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Diagrams shown in fig. [I]p 



27e 27(1 - s) 2 
2s 



7- 15s + 8s 2 -nis(l -3s) 



27(1 - s) 3 
iP(s) 4 



(l - 4s + 3s 2 - 27R*(1 - 2s)) ln(s) + (1 - 2s) ln 2 (s) 



(52) 



2 / \ 2 1 
r 4 ln(s) - 5 - vr 2 H — 

27e 2 27e V V ; / 9 27(1 - s) 2 



7- 10s + 3s z 



-4?ri(6 - 13s + 5s 2 



4 



ln(s) 



6 - 19s + 18s 2 - 5s 3 



27(1 - s) 3 

+2?ri(l - 4s + 6s 2 - 2s 3 ) - (1 - 4s + 6s 2 - 2s 3 ) ln(s 
Diagrams shown in fig. [I]p 

<P(s) = J- - |[ ( 12 ln ( S ) - 49 + 48 C(3) " 127TZ 
Diagrams shown in fig. [^1 

yP(s) = — + — ( 37 + 12m) + ^— ln(s) 

v ; 243e 729 V J 243(1 - s) K ' 

4 2 / , _ \ 1 

■ 121n(s) - 19 - 127ri 



(53) 



(54) 



243e 2 729e V 



2187 



(55) 

72 Li 2 (s) + 72 ln(l-s)ln(s) 



+36 ln 2 (s) + 463 - 90vr 2 - 12ttz(12 ln(s) - 19) 
Diagrams shown in fig. 

Mis)- 



729(1 - s) 



ln(s)(13- 19s) (56) 



1 ° + L ( 10 - 27s + 33s 2 - 10s 3 )arccot(Vz- 1) 



243e 243(1 - sf 



107 + 452s - 696s 2 + 464s a - 113s 4 - 6vr z s 



729(1 - s) 4 

+216s arccot 2 (V^ - 1) 



(57) 



4 2 
+ 



243e 2 729se 

2v^T(2 + s 
729~i 



(24y/z- 1(2 + s)arccot(Vz - 1) - 48 - 19s 



6tt f ln(xi) + ln(x 2 ) - 2 ln(z - 1) 
+i(24Li 2 (x 2 /x!) - 12Li 2 (-a; 2 /x 1 ) - 15 ln 2 ^,/^) + 121n 2 (x 2 ) 
-12 ln(xi) \n(x 2 ) + 24 ln(xi) ln(x 4 ) - 24 ln(x 2 ) ln(x 4 ) - 5vr' 
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^6 + s — s 2 ) ^47rarccot(V-2 ; — 1) — ln 2 (x 2 /a;i) j 



81s(l - s) 

- — ( 1056 - 1541s + 130s 2 + 463s 3 ) 

1 — s) V / 



2187s(l - s) 2 



(32 - 65s - 6s 2 + 40s 3 - 19s 4 )arccot( v / I 3 T) 
6^3 - 8s + 2s 2 )arccot 2 (v / ^ - 1) 



729s(l - s) 3 
4 

+ 729s(l - s) 4 

-7T 2 (^27 - 67s + 28s 2 + 36s 3 - 34s 4 + 7s 5 )] 

C One— loop matrix elements to O(e) 

The one-loop matrix elements 

(d£+r ((^one-loop = ^; l0 ° P(7) (a)t r ee + Wtree 

are needed to 0(e) for the operators ©1,2,4,11, 12- The results are 

p l-loop(9) _ 3 i_i oop (9) 
2,m — ^ ^l,u 

Fl^' 7) = (f ) { " 5 + f (V^TarccottV^T) - l) } 

F i-ioo P (9) _ f ~ £ r i6 



11 \ fi 2 J { 27e 

8 r 

t 



81s 

£(2 + s)y^^T 

8li 



6(2 + s)y/z - larccot( v / ^ - 1) + 3sln(s) - 12 - 4s - 3vris 

1 (48 U 2 (x 2 /x l ) + 24 ln 2 ^) - 12 ln 2 (x 2 ) + 12 In 
+24 ln(xi) Info) - 48 Info) Info) - 3 ln 2 (^ - 1) - 5vr 2 ) 
+6tt ^6 Info) + 2 Info) + 2 Info) + 2 Info) - hi(z - 1)) 



32e(5 + s) v / ^T , , -. 

H arccot(V-2 — 1) 

Olo 



As 
243i 



120 + (52 + 12m - 9tt 2 )s - 6(2 + 3™)s ln(s) + 9s ln 2 (s) 
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p l-loop(9) _ 3 i_i oop (9) , , 
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